Entrainment coefficient and effective mass 
for conduction neutrons in neutron star 
crust : simple microscopic models 

Brandon Carter", Nicolas Chamel", Pawel HaenseP'^ 

On 
O 

o 

CN ! Observatoire de Paris, 92195 Meudon, France 

^ '■ (Brandon.Carter@obspm.fr, Nicolas.Chamel@obspm.fr) 

^ N. Copernicus Astronomical Center, Warsaw, Poland 
(N| ; (liaensel@camk.edu.pl) 



X 



February 8, 2008 



> 

in 
o 

^ ■ Abstract. In the inner crust of a neutron star, at densities above the "drip" 

Q . threshold, unbound "conduction" neutrons can move freely past through the ionic lat- 
tice formed by the nuclei. The relative current density 77.* = 777)* of such conduction 
neutrons will be related to the corresponding mean particle momentum pi by a propor- 
O I tionality relation of the form 77* = JCp^ in terms of a physically well defined mobility 
^ ■ coefficient JC whose value in this context has not been calculated before. Using methods 
from ordinary solid state and nuclear physics, a simple quantum mechanical treatment 
based on the independent particle approximation, is used here to formulate JC as the 
phase space integral of the relevant group velocity over the neutron Fermi surface. The 
I result can be described as an "entrainment" that changes the ordinary neutron mass 777 
to a macroscopic effective mass per neutron that will be given - subject to adoption of a 
convention specifying the precise number density 77 of the neutrons that are considered 
to be "free" - by 777^ = n/JC . The numerical evaluation of the mobility coefficient is 
carried out for nuclear configurations of the "lasagna" and "spaghetti" type that may 
be relevant at the base of the crust. Extrapolation to the middle layers of the inner 
crust leads to the unexpected prediction that 777^ will become very large compared 
with 777 . 
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1 Introduction 



The main purpose of this article is to show how a mean field treatment of neutron star 
crust matter can be used to address the previously unsolved problem of evaluating a 
quantity, namely the relevant neutron mobility coefficient KL , that is essential for the 
astrophysical applications that will be described in a separate article pp. In terms of the 
relevant number density n of effectively unbound neutrons, this coefficient determines 
a corresponding effective mass 772^ = n/K, that characterises their average motion 
on a macroscopic scale (meaning one that is large compared with the spacing between 
nuclei) and that we therefore refer to as the macro mass, to distinguish it from the 
microscopic effective mass, 772® say, characterising the dynamics of the neutrons on 
subnuclear scales. Whereas Tfl® is well known to be typically rather smaller than the 
ordinary neutron mass Tfl , we reach the previously unexpected conclusion that there 
is likely to be a strong "entrainment" effect whereby the macro mass 772^ will typically 
become large, and in some layers extremely large, compared with m . 

A secondary purpose of this article is to draw the attention of nuclear theorists to 
the potentialities of the almost entirely unexplored branch of theoretical astrophysical 
nuclear physics that needs to be developed for this and many other purposes. The 
only relevant work of which we are aware so far is that of Oyamatsu and Yamada |2] , 
who appear to be the only ones to have taken proper account of the neutron scattering 
by the nuclei inside the inner crust by the use of appropriate Bloch type periodicity 
conditions of the kind commonly employed for the treatment of electrons in ordinary 
terrestrial solid state physics. Their treatment however was restricted to a simple one 
dimensional model. 

Under the conditions of ordinary terrestrial solid state physics, and even at the much 
higher densities characterising the matter in a white dwarf star, long range electric 
forces keep the nuclei so far apart that, in so far as the much stronger but short range 
nuclear interactions are concerned, each individual nucleus can be treated separately 
as if it were isolated. Until quite recently [2], such a separate treatment of individual 
nuclei (considered as if isolated each in its own cell - with Wigner Seitz type boundary 
conditions) has been used in nearly all quantum mechanical calculations on neutron 
star crust matter since the pioneer work of Negele and Vautherin That kind of 
approximation is fully justifiable in the outer crust , where the densities are not too much 
greater than those found in a white dwarf. However such a treatment can no longer be 
considered entirely satisfactory in the inner crust, meaning the part with density above 
the "neutron drip" threshold at about 10^^ g/cm^, where there are unconfined neutrons 
that travel between neighbouring nuclei, which thereby cease to be effectively isolated 
from one another. 

While desirable for accuracy throughout the inner crust, a proper collective rather 
than individual treatment of the nuclei becomes not just desirable but absolutely essen- 
tial for treating the problem with which Oyamatsu and Yamada |2| jl] were concerned, 
namely that of the nuclear matter inside neutron star crust. Such a treatment is also 
essential for treating the problem with which the present work is concerned, namely 
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that of stationary but non static configurations in which a neutron current fiows rela- 
tive to the lattice formed by the nuclei, something that obviously can not be discussed 
in the usual approach that treats the nuclei as if they were isolated in individual (e.g. 
Wigner Seitz type) boxes. 

The flow of neutrons is treated here as a perturbation of a zero temperature ground 
state characterised just by the location of the relevant Fermi surface in momentum 
space. We thereby obtain provisional rough estimates of the relevant mobility coefficient 
which suggest that (unlike what occurs in the fiuid core and on a microscopic scale) 
the macro mass that effectively characterises the neutron motion on a macroscopic 
scale can become very large compared with the ordinary neutron mass m , particularly 
in the middle part of the conducting layer, for which three dimensional numerical results 
will be presented in a follow up article jB]. The present article deals more specifically 
with simplified rod and plate type models that are relevant near the crust core interface, 
where the mass enhancement will be less extreme. 

Bulgac, Magierski and Heenen have recently pointed out the importance of shell 
effects induced by unbound neutrons in neutron star crust by evaluating the Casimir 
energy for neutron matter in the presence of inhomogeneities from a semiclassical ap- 
proach and more recently by performing a Skyrme Hartree-Fock calculation with or- 
dinary periodic boundary conditions (see [H |Hj and references therein). However this 
kind of boundary conditions does not properly account for Bragg scattering of dripped 
neutrons and is only a particular case of the more general Bloch type boundary condi- 
tions. 

In the absence of any previous quantum mechanical calculation whereby nuclei on 
a crystal lattice are treated collectively (apart from the ID calculation previously men- 
tionned j2]), even at the simplest level of approximation, we shall adopt the simple 
model suggested by Oyamatsu and Yamada, supplemented with Bloch type boundary 
conditions, in order to estimate the effective neutron mass . This model treats the 
neutrons as independent fermions subject to an effective potential. We wish to draw the 
attention of nuclear theorists to the problem of including Bloch boundary conditions in 
more sophisticated approximation schemes as a challenge for future work. 

In the mean time, experience with the analogous problem of electron transport 
in ordinary solid state physics suggests that the results obtained from the Oyamatsu 
Yamada type treatment used here should not be too bad as a first approximation. 
Further encouragement comes from our own recent attempt to take up the challenge 
of allowing for coupling by an appropriate adaptation of the standard BCS pairing 
theory on which the prediction of neutron superfluidity (in the relevant low to moderate 
temperature range) is based: the upshot [0] is that (although it is essential for the 
inhibition of resistivity) as far as the "entrainment" phenomenon is concerned the effect 
of the ensuing pairing "gap" will not be very large nor very difficult to calculate. 
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2 Microscopic description of conduction neutrons 
in the inner crust 



2.1 Single-particle Schrodinger equation 

The basic principle of the conductivity model we wish to adapt from ordinary solid 
state theory to the context of a neutron star crust is that the "conduction band", 
and perhaps also some of the highest "confined" levels, can be analysed within the 
independent particle approximation, in terms of energy eigenstates for a single particle 
described by a Bloch type wave function Lp , satisfying the Bloch periodic boundary 
conditions as discussed below. In the rest frame of the crust, with respect to which 
the system will be assumed to be in stationary equilibrium, the single particle wave 
function will be taken to be governed by a Hamiltonian operator 1~L that is given by 

n = -n^\/^^Vj + v , (2.1) 

2m® 

where 7*-^ is just the space metric, while V is the single particle potential, and TTl® 
is the relevant local effective mass parameter. 

The potential V and the effective mass TJl® will be periodic in the case of a regular 
crystalline solid lattice, though not for a fluid configuration such as will be relevant 
at higher temperatures. The value of the effective potential V (and the associated 
deviation of the effective mass TJl® from m ) is supposed to allow not just for the 
attraction of the confined protons in the nuclei but also for the mean effect of the other 
fermions, which are electrons in the familiar solid state applications, but neutrons in 
the context under consideration here. 

2.2 Boundary conditions 

The preceding considerations apply even to disordered (glass like or liquid) configura- 
tions, but in order to proceede we shall now restrict our attention to cases for which 
the nuclei are assumed to be fixed and locally distributed in a regular crystal lattice 
configuration, in which an elementary cell consists of a parallelopiped of volume Vcell 
say spanned by a triad of basis vectors e^ labelled by an index with values a = 1, 
2, 3, so that the system is invariant with respect to translations generated by vectors 
of the form 

T = re« (2.2) 

for integer values of the coefficients (in the following summation over repeated 
indices is assumed). This so called adiabatic or Born-Oppenheimer approximation is 
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justified by the large difference between the neutron and nucleus masses, since typically 
each nucleus contains several hundred nucleons |lUj . 

It follows directly from the well-known Floquet-Bloch theorem ^I] that the single 

particle wave function has to satisfy the following boundary conditions 

(^k{r + T} = e^'^-V{r}. (2.3) 

where k is the Bloch momentum covector. 

It will therefore suffice to solve the Schrodinger equation just inside a single elemen- 
tary cell. Instead of using a primitive cell of parallelopiped, it is convenient for many 
purposes to work instead with the Wigner-Seitz (W-S) cell defined as the set of points 
that are closer to a given lattice node than to any other. Such a cell exhibits the full 
symmetry of the lattice. Its shape is determined by the crystal structure, a polyhedron 
in general, for instance a cube in a simple cubic lattice. This exactly defined W-S cell 
should not be confused with the widely employed eponymous "W-S approximation" 
which consists in replacing this cell by a sphere (or more generally any convenient 
cell that simplifies the analysis). 

The Bloch momentum k takes values inside the first Brillouin zone (B-Z), which 

is the W-S cell of the reciprocal lattice whose nodes are located at K = ia^"" for 

integer values of the coefficients ia where the dual basis vectors l'^ are defined by the 

following scalar products 

r-e, = 27r(5%, (2.4) 

the 27T normalization factor being introduced for convenience. 

It must be emphasized that the single particle wave function will have to satisfy the 
relation ()2.Hj) between two opposite faces of the cell. This means in particular that the 
Schrodinger equation has to be solved for each wave vector k inside the first B-Z. The 
ordinary periodic boundary conditions on the cell, which have been recently applied by 
Magierski and Heenen [S], are thus only a restricted subset of solutions, namely those 
with k = 0. 

For each momentum k , there exists only a discrete set of energy eigenvalues satis- 
fying the Bloch boundary conditions. The single particle energy spectrum is therefore 
a collection of sheets in momentum space, each sheet usually being referred as a band 
(labelled by the principal quantum number). 

3 Microscopic dynamics in mean field of lattice 

3.1 From microscopic to macroscopic observables 

The linearisation involved in the neglect of direct two body or many body interactions 
in such a mean field treatment excludes allowance for higher order effects such as the 
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pairing responsible for a superfluid energy gap, but as discussed in a separate article 
jH], this limitation should not matter too much for the evaluation of the basic equation 
of state of the fluid, since the main effect of superfluidity is not to modify the equation 
of motion but just to restrain the class of the admissible solutions (by allowing only 
those that are irrotational). 

The main limitation on the use of such a linearisation is that it makes sense only 
for configurations that do not differ too much from the static reference configuration 
on which the estimation of the effective potential energy function V is based. 

3.2 Fermi surface in ground state configuration 

The zero temperature configuration is obtained by minimising an energy density U 
subject to the constraint of a given value of the neutron density TT-n 5 defined by 

"""^Jj^,- (3.1) 

The energy density is expressible in terms of the single particle energies (in this work, 
we shall use braces instead of ordinary brackets for functional dependence, in order to 
avoid possibly confusion with simple multiplication) 

nipi, = S{k}ip^, (3.2) 

by an expression of the form 

U = U{0} + l£{k}—^, (3.3) 

where, as a standard postulate, the contribution U {0} is ignored in the minimisation 
procedure. This means that all single particle states of energy up to but not beyond 
some particular Fermi level, ^ say, are occupied. This Fermi energy specifies the 
"Fermi surface", namely the locus, say, in momentum space, where 

£{k} = fi. (3.4) 



It must be emphasised that the Fermi surface will in general consist of disconnected 
pieces unlike the homogenous case in which the Fermi surface is simply a sphere. From 
the property that ^{k} = — k}, the ground state is thus completely symmetric 
in the sense that its total momentum density / d'^k /c* vanishes. This is a generalized 
version of Unsold theorem according to which the all electron wave function of a 
closed shell atom is spherically symmetric. 
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3.3 Minimal conducting configurations 

Our main concern in the present analysis is with current carrying configurations, as 
characterised by some given value of the physically well defined number current density 
with components given by 



in which the transport velocity V is given by the formula 

The kind of current carrying configurations that will presumably be relevant in the 
low temperature limit will be those for which the energy density U is minimised for a 
given value of the neutron density subject to the further constraint that the number 
current density also has the prescribed value 77-*. It is evident that this constrained 
minimisation condition requires that the occupied states should consist just of those 
subject to an inequality of the form 

£<|1+P^V\ (3.7) 

for some fixed covector whose constant components Pi have been introduced as La- 
grange multipliers. The phase space volume specified in this way will have a boundary 
characterised by the condition 

S = ii + p,v\ (3.8) 

which specifies a modified Fermi surface, S say. The constraint of having a finite 
current breaks the symmetry and will thus necessarily lead to a new ground state with 
a non vanishing total momentum density vector / d'^k fc* 7^ . 

3.4 The uniform displacement of the Fermi surface 

To justify the use of a Schrodinger type Hamiltonian Ti, involving an effective potential 
of the same form as in the zero current case, we need to assume that (as will presumably 
be the case in the relevant context of pulsar glitches) the current density 77* is small 
enough to be treated by linear perturbation theory. This means that the Lagrange 
multiplier Pi should itself be considered just as a first order perturbation, having zero 
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Figure 1: Sketch of energy against wave number, showing uniform displacement at- 
tributable to current. 

value Pi = in the static unperturbed configuration, and Pi = 6pi in the perturbed 
current carrying configuration. For a given (unchanged) value of the chemical potential 
p , the difference between the new value ()3.8p and the old value ()3.4|) of the Fermi energy 
level S will be given to first order by = 6[pi V^) = 6pi and thus in terms of 
the perturbed value of pi simply by 

6S=p^v\ (3.9) 

This change in energy can be interpreted as being attributable to a phase space dis- 
placement 6ki. Since this change will be given, according to ()3.6|) , by 

6S = nv'6k,, (3.10) 

it can be seen that the simplest possibility is to take the displacement to have the 
uniform value given just by 

hSki=p^, (3.11) 

as illustrated on figure [U 

3.5 Relation between current and momentum 

The relation (jH.llll means that (in the infinitesimal limit) the Lagrange multipliers Pi 
can be physically interpreted simply as components of a uniform pseudo momentum 
displacement of the occupied phase space region, and of its "Fermi surface" boundary. 

The effect of this displacement on the Fermi surface element dS^^ as given - in 

terms of the corresponding surface measure element, dS^ - by 

dS'/ = {vyv)dS^, (3.12) 
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will be to sweep out an infinitesimal phase space volume element given by 

hd^k = p,dSl. (3.13) 

It follows that, for the integral over the occupied region of any phase space function 
/ , the difference between the value for the displaced (conducting) configuration and 
the value for the non conducting reference configuration will be given (to linear order) 
by the formula 

h6jfd'k = p,fjdSl, (3.14) 

in which it is to be understood that the integral on the right is taken over the entire 
Fermi surface. 

We have already pointed out that the energy function S will be symmetric with 
respect to the origin, hence it follows that the Fermi surface element dS^ will be 

antisymmetric about the origin so that its unweighted integral dS^ will self cancel 
to give zero. We can thus see from ()3.14j) that for any phase space function / that 
is symmetric about the origin the corresponding integral will be unaffected by the 
displacement, i.e. we shall get Sjfd^k = . This contrasts with the case of an 
antisymmetric function, for which it is the unperturbed integral that will vanish, i.e. we 
shall have J f d'^k = in the static reference configuration, so that the corresponding 
value for the conducting configuration will be given by / / d^k = 6 J f d'^k . 

The antisymmetric case is illustrated by the application with which this work is 
principally concerned, namely the current density 71* given by ()3.5j) . Since we have 
= in the reference configuration characterised by J9j = , we shall have = dn^ 
in the conducting configuration characterised by Pi = 6pi . Thus by substituting 
for / in (nmi we see that in the linearised limit the current will be related to the 
pseudo momentum displacement Pi by a relation of the form 

n' = ]C^p^, (3.15) 

in which the symmetric tensor A^*-' is given as an integral over the Fermi surface by 

1 r v'^V^ 

/C*^ = — dS^. (3.16) 

Before continuing, it is to be remarked that this tensor lO'-^ is interpretable as 
being proportional to the zero temperature limit of the electric conductivity tensor d*-^ , 
defined in the usual way by J* = cr^'^Ej where Ej is the relevant electric field and 

J* the electric current density, by a relation of the form cr*-^ = re^AT*-^, where T is 
the relevant relaxation time and e the electric charge per particle. 
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In a entirely disordered (liquid or glass like) state the form of the energy function 
8 and the consequent location in phase space of the Fermi surface will be difficult to 
evaluate theoretically, but there will be the partially compensating simplification that 
the result will automatically be isotropic. In the mathematically simpler case of a cubic 
crystalline lattice that is expected ^3] to occur in a neutron star crust, this tensor KL^^ 
will have the isotropic form 

r' = jcY^, (3.17) 

in which the scalar coefficient will evidently be given by 

}C = l-f,j}C'r (3.18) 

However that may be, it can be seen from that the scalar coefficient JC will be 

given by a Fermi surface integral of the simple form 

In terms of this integral, the relation between current and momentum will be given 

by 

= JCY^pj . (3.20) 



3.6 The (cut-ofF dependent) concept of the effective mean mass 

In order to relate the formula ()3.19|) to expressions used elsewhere in the literature, it is 
to be noted that we can introduce a mean velocity vector and a conduction neutron 
density n such that the current density can be expressed by 

n' = nv' . (3.21) 



Subject to the isotropy condition ()3.17|) , the (pseudo) momentum covector pi will 
be expressible in terms of the corresponding covector = ^ijV-^ in the simple form 

Pi = m^Vi , (3.22) 

in which the effective mean particle mass - which we refer to as the "macro mass" 
in order to distinguish it from the locally effective '"micro mass" Ul® - is defined in 
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terms of the integral (j3.19p by the relation 

1 



(3.23) 



n 



It is however to be remarked that whereas the specification of quantities such as Pi 
and JC is physically unambiguous, the specification of the effective mass , like 
that of the mean transport velocity , depends on how many "conduction states" are 
counted in the definition of the number density n , and how many are left aside as 
dynamically inert "confined" states. 

In the application to a neutron star crust, the situation is somewhat simpler than 
is usual in ordinary solid state physics because the relevant single particle Hamiltonian 
fl2.1|) will usually involve an effective potential function V that tends rapidly [T3] 
towards an almost exactly uniform value outside the ionic nuclei to which all the protons 
and some of the neutrons are effectively confined. It should be noticed that the single 
particle wave function of bound states will be vanishingly small at the W-S cell boundary 
and therefore those states will not be sensitive to the Bloch phase shift. The resulting 
single particle energies will thus be nearly independent of the Bloch momentum which 
means that the associated group velocities will be very small. The single particle energy 
spectrum can therefore be decomposed into a subset consisting of such confined states, 
and a remainder that will realistically be describable as "conduction" states. The 
separation between those two different subsets is not entirely sharp (borderline states 
are commonly refered to as valence states in ordinary solid state physics) and one has 
to rely on a more or less arbitrary convention. 

In most layers of a neutron star crust, the most convenient possibility is usually to 
take the uniform value of V outside nuclei, which can be taken as the energy origin, to 
specify a corresponding energy range 8 > characterising "conduction" states. This 
specification may however be ambiguous in the bottom layers of the inner crust where 
nuclei are very close to each other. In order to deal with such cases, we shall adopt 
the convention, which agrees with the previous one when nuclei are very far apart, that 
conduction states are defined as states whose energy is larger than the maximum value 
of the potential. 

It is to be remarked that in the uniform limit whereby the periodic potential V 
and the local effective mass m® are simply constants, the resulting Fermi surface is 
a sphere and the evaluation of the mobility scalar is straightforward. In this case, the 
associated macro mass is just equal to the micro mass = m® provided all neutron 
states are counted as conduction states. 

4 Quantitative estimates of the mobility tensor 
4.1 Bottom layer of the inner crust 

We shall focus for simplicity on the bottom of the inner crust near saturation density 
of the order of 10^^ g/cm^ where the transition to homogeneous fluid neutron matter 
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spaghetti 



lasagna 



Figure 2: Nuclear configurations for liexagonal spaglietti (rod) and lasagna (slab) lat- 
tices. Black - nuclear matter; grey - neutron gas. From with kind permission of the 
author. 

takes place, to evaluate the mobility scalar and the effective mass. In this region the 
crystal lattice is not expected to substantially alter the transport properties such as 
the mobility tensor. The aim of the following sections is not to provide astrophysically 
important information but rather to give some insight that may be valuable in the more 
general cases that will be dealt with in future work. Near the base of the crust, the 
nuclei are so strongly deformed by their neighbours that they may adopt non spherical 
shapes ^H], sometimes idealised by ID or 2D configurations such as slabs and rods 
respectively, which greatly simplifies the analysis. These "exotic" crust phases are 
illustrated on figure El taken from the work of Oyamatsu j3] . 

In views of the exploratory character of the present work, we shall use a simple model 
for the single particle equation. In particular, we shall drop the condition of strict self- 
consistency, and shall use a mean field model based on the work of Oyamatsu jH E] . 
He calculated the structure of the ground state of the inner neutron star crust using a 
phenomenological energy-density functional, fitted on the first hand to the smoothed 
nuclear masses and charged radii of laboratory nuclei on the j3 stability line and on the 
other hand adjusted on the equation of state of pure neutron matter from the variational 
calculations of Friedman and Pandharipande [3] using two body as well as three body 
nucleon-nucleon interactions. He further investigated with Yamada the importance of 
shell effects with a single particle Schrodinger equation in the W-S approximation [2] 
(the reader is refered to those two papers for the details). 

Unlike the calculations carried out by these authors, we shall apply the ID or 2D 

Bloch type boundary conditions and we shall not make any approximations about the 

shape of the W-S cells which partition the crystal. The 3D configurations require specific 

numerical techniques and will be discussed in a separate work (Hj. For consistency, the 

lattice parameters (whose values are found in ^) will be defined so that the volume 

of the exact W-S cell is equal to the volume of the approximate cell. Within a W-S 
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cell, we approximate the single particle potential V by the potential U of Oyamatsu 
and Yamada Pj. A zeroth order approximation Uq of the potential is obtained by 
differentiation of the potential part of an energy density functional V^Tlj^, rip} so that 

Specifically we shall use the parameter set of model I of this paper. Then, the effect 
of the finite range of the nucleon-nucleon interaction is taken into account by using 
the folding of the potential U q with a gaussian smearing function of width t\. This 
potential is supplemented with a spin-orbit coupling term, which is parametrized as a 
functional of the density gradients. The gaussian width and the parameters of the spin- 
orbit potential are adjusted so as to reproduce the correct sequence of single particle 
energy levels of ^°®Pb. The potential U thus obtained varies only near the nuclear 
surface. Therefore, close to the W-S cell boundary U is constant. This enables us to 
express the periodic potential V acting on a neutron moving in an infinite lattice of 
nuclei as 

V{v] = Y.U{t-T] , (4.2) 

T 

where the sum goes over all nuclei occupying the lattice sites, i.e. T = l^'G.a , where 
are integers. The potential V^{r} thus possesses all the symmetries of the crystal 
lattice. We neglect small changes in U due to passing from a single W-S cell to an 
infinite lattice. The energy origin is taken as the value of the potential V outside 
nuclei. Following Oyamatsu and Yamada [2] we do not include momentum dependent 
terms in the single particle potential. Consequently the microscopic effective neutron 
mass in this model is just equal to the bare neutron mass Ul® = m . We have ignored 
the spin-orbit coupling since it is found to be about one order of magnitude smaller 
than the central potential [2]. 

The single-particle Schrodinger equation is solved here by the Rayleigh-Ritz varia- 
tional approach whereby the expectation value of the Hamiltonian is minimized subject 
to a normalization condition. The single particle wave function is expanded into plane 
waves defined by 

^k{r} = E ^k{K} , (4.3) 

V ^^cell K 



with the normalization 



El^k(K)P = l. (4.4) 

K 
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This expansion into plane waves is actually exact, it merely makes more apparent 
the periodicity. The approximation lies in the fact that the summation needs to be 
truncated for practical calculations at some cut off energy ^cutoff (which thus makes 
this method adequate only for slowly varying potentials as in the present case), i.e. so 
as to include all reciprocal lattice vectors satisfying 



2m 



< ^cutoff • (4.5) 



The group velocity is found according to the Hellmann-Feynman theorem ^7] to be 
equal to 

v = — + E^^k K 4.6 
m K 

This last equation shows that deviations from the homogeneous case arise from the 
spreading of the Bloch wave packet ()4.3|) . 

For each total neutron density , the Fermi energy /i is determined as an integral 

over the first B-Z by 

^n = 7T^E/ d'k^{f,-S^}, (4.7) 

where we have introduced the Heaviside unit step distribution defined by f^jx} = 1 
if X > and zero otherwise (the factor of 2 is to account for the spin degeneracy). 

The conduction neutron density is defined by the occupied single particle states whose 

energy is positive, i.e. 8^ > , namely 

2 



o x3E/ d^fc^A^-^jn^J. (4.8) 



n 

a "'BZ 



It will be instructive to compare the Fermi surface area with that of a non interacting 
neutron gas of density 77.^ , which is given by 

^gas = 47r(37r2nn)'/' . (4.9) 



4.2 "Lasagna" phase 

In this model, the slab shaped nuclei or "lasagna" are parallel to each other and equally 
spaced by a distance a . Hence the single particle potential V is periodic along one 
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direction only, say the Z axis, and takes a constant value in the other two dimensions. 
The neutron band structure with a square well potential (whose analytic solution was 
found a long time ago by Kronig and Penney [IHl EH] ) was the only case discussed by 
Oyamatsu and Yamada [2j. 

The single particle wavefunction can be factored in the form 

^k{r} = 0t.{4e^^'-"+'^^\ (4.10) 

in which the reduced wave funtion is the solution of a one dimensional equation 
of the form 

^''^ ' v{z}<\>,Xz} = E{kMkXA. (4.ii: 



2m dz'^ 
satisfying the boundary conditions 

0,,{z + a} = e^^->,J4. (4.12) 



The first B-Z is defined by the set of Bloch wave vectors such that — 7r/a < kz ^ 
7r/a. The B-Z in this peculiar case has an infinite extent as a result of the continuous 
translational invariance along planes parallel to the "lasagna". The single particle 
energy is thus given by an expression of the form 

i:„{k} = 5^^^|±^ + aM, (4.13) 

where Ot is a band index. It follows that the neutron group velocity components 
parallel to the slabs coincide with the group velocity of a non interacting neutron gas 
while only the component perpendicular to the slabs is affected. It can be shown from 
group theoretical arguments, that this component of the group velocity must vanish at 
the Brillouin zone edge 120)5 i-e. at = :^Tl / a . More generally whenever the Bloch 
wave vector lies on a symmetry plane, the group velocity component normal to the plane 
vanishes . It can also be shown that the energy bands do not cross |2Uj (nevertheless 
bands may touch "accidentally" for some specific choice of potential, for instance a 
constant one). 

The total neutron density is given by 



m „ /-Tr/a 



TT 



E r dA:,(^ - Eo)^{^i - £ J . (4.14) 



The mobility along the "lasagna" planes coincides with that of the uniform neutron 
gas with the same total neutron density: 

k}\ = )c^^ = )cyy = ^-^. (4.15) 

m 
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The mobility in directions perpendicular to the "lasagna" may however differ from 
that of the non interacting neutron gas: 

(4.16) 

Since the "lasagna" configuration is strongly anisotropic, it is more appropriate to 
define a transverse effective macro mass by 

mi ^ ^. (4.17) 

The conduction neutron density is expressible as 

n = nn + 777^ E L ' TT^^^i?!^ - Eo^-di-Eo^^k^ . (4.18) 



4.3 "Spaghetti" phase 

In the spaghetti model, the crystal is composed of cylinder shaped nuclei arranged on a 
two dimensional lattice (the nuclear "spaghetti" are assumed to be parallel to each other 
along say the z axis). Since the single particle potential around one isolated nucleus 
only depends on the distance to the rod like nucleus, the corresponding contribution to 
the crystal potential does not depend on z hence the wave function can be factored as 

¥^k{r} = 0,,,,,{x,^}ei^^^^ (4.19) 

The two dimensional wave function obeys a Schrodinger equation of the form 

(4.20) 

The energy is thus decomposible as 

^«{k} = e^{k,, ky} + ^. (4.21) 
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The total neutron density is given by 



4 



/^/ ■'R7. 



in which the factor of two arises from the restriction that fc^ > and where the 
integration is carried out over the 2D first B-Z illustrated on figure lU 

It is readily verified that the mobility tensor component along the cylindrical nuclei 

is merely equal to the mobility component of the non interacting neutron gas: 

/Cll =}C'' = —. (4.23) 
m 

The other components perpendicular to the "spaghetti" depend on the neutron-crystal 
interaction. It is convenient to defined a mean transverse mobility by 

fm 2 r dk^dky f /dSgY {d£a\A .r _ 

V 2 {27,fn^^L [ydkj ^ \dky) ) 

(4.24) 

since it involves the integral of a completely symmetric function that can thus be 
factored into one irreducible domain of the first B-Z fU] . 

The Fermi surface area, the mobility tensor and the density involve integrations over 
the two dimensional first B-Z of functions weighted by a Fermi distribution. We have 
followed an idea due to Gilat and Raubenheimer for three dimensional crystals [22]; 
and translated it into the two dimensional case. First of all the zone is decomposed 
into small identical cells within which the single particle energy is linearly extrapolated 
from the value of the energy and its gradient at the center of the cell. The integration is 
then performed analytically inside the cell. The integral is approximated by summing 
the contribution from each cell (properly weighted whenever the cells overfill the zone). 
Unlike the case of band crossing, extrapolation fails for band "kissing" where the energy 
gradient varies rapidly in the cell (see the schematic pictures on Figure O)). This problem 
becomes more and more acute as the number of bands to be included in the integration 
is increased. However those induced systematic errors will tend to vanish as the number 
of microcells is increased. Whereas the integrand could be integrated analytically inside 
each microcell, we found that it is better to take it as constant. The reason again lies 
in the band "kissing" problem. The integrand typically depends on the momentum 
via the single particle energy which is extrapolated. Near a band kissing region, the 
extrapolation errors in the integrand are integrated which lead to unstable results. 

We have considered two lattice types: square and hexagonal crystals (whose recip- 
rocal lattices are also square and hexagonal respectively). The point group 'P (the 
set of point symmetries which send the lattice into itself) of the hexagonal lattice is 
Cgy (Schonflies notation |23|) whose order is equal to = 12, thereby reducing 
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band crossing bandltissing 

Figure 3: Sketch of energy against wave number for crossing and "kissing" scenarios. 

integrations of completely symmetric functions over the entire two dimensional first 
B-Z to 1/12*^ of the zone. We mention that for this structure the lattice spacing a 
(shortest distance between any two lattice points) is not equal to the value given by 
Oyamatsu a^^ but is given by 



a 



\ 



(4.25) 



Likewise, the first B-Z of a square lattice whose point group is C^y (lattice spacing 
a = ) can be partitionned into 8 irreducible domains. These considerations are 
illustrated in Figure |3] with the conventional labelling (first B-Z in red, one irreducible 
domain in green). The most adequate cell for integrations is thus a rectangle since the 
irreducible B-Zs are rectangular triangles. The partition of this domain into microcells 
is thereby straightforward. 

These two groups, Cq^ and C4^y , contains two dimensional irreducible representa- 
tions 1^ which means that unlike the ID case, energy bands may cross each other (20]. 
The neutron band structure along high symmetry lines is shown on Figures El and IHl for 
energies in the vicinity of the Fermi energy and for the lattice spacing a^^ = 27.17 fm . 



The conduction neutron density is equal to 



n = n 



4 f 

n - ^27i)3n ^ Jbz V-2mec,^{A^ - £a}^{-£a}dk^dky . (4.26) 



The transverse effective mass is defined by 



n 



(4.27) 
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Figure 4: First Brillouin zones and irreducible domains for square and hexagonal vari- 
ants of spaghetti (rod) lattice. 
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Figure 5: Neutron band structure around the Fermi energy ^ , along directions shown 
on figure m for square variant of spaghetti (rod) lattice (Oq = 27.17 fm). 
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Figure 6: Neutron band structure around the Fermi energy /i , along directions shown 
on figure m for hexagonal variant of spaghetti (rod) lattice with a^^ = 27.17 fm . 
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Figure 7: Convergence of computed effective mass Tflj^ /m as a function of cell number 
for an hexagonal lattice with a^^ = 27.17 fm . 
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Figure 8: Effective mass of conduction neutrons 772^ /m as a function of lattice spacing 
{(^Qy ) below and above phase transition from lasagna (slab) regime to spaghetti (rod) 
regime. 



The convergence of the integration scheme based on a decomposition into rectan- 
gular cells is illustrated on figure [7| for the hexagonal lattice with the lattice spacing 
a^^ = 27.17 fm . The convergence is much faster for the density than for the other 
quantities due to the absence of the singular square root integrand. 



5 Discussion 



The results concerning the macroscopic effective neutron mass are illustrated on figure |H1 
and numerical values can be found in the appendix. The macroscopic effective neutron 
mass appears to be increased compared to the ordinary neutron mass. The figure 
also shows that this renormalisation of the neutron mass is mostly significant at low 
densities and becomes negligible at higher densities where nuclei nearly merge into a 
uniform mixture. The deviations of the effective neutron mass from the bare one can be 
understood in terms of modifications of the Fermi surface from a sphere. In particular, 
as a result of the opening of band gaps the Fermi surface is not smooth but contains 
many holes as schematically illustrated on figure IHl Since the enclosed Fermi volume 
only depends on the density, this means that the Fermi surface area for a given density 
is reduced as compared to the corresponding sphere as shown on figure ^1 (see the 
appendix for the numerical values). In the case for which the Fermi volume is equal to 
the volume of the first Brillouin zone, the Fermi sphere can be distorted and torn such 
that the resulting surface area simply vanishes while the volume remains finite. Such a 
situation occurs in ordinary electric insulators. 

In the present work we have neglected the spin-orbit coupling which is a small per- 
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Figure 9: Example of energy contours in extended Brillouin zone for square lattice. For 
sufficiently small Fermi energies the contours are approximately circular as in the free 
particle case, while at higher energies, the contours consist of disconnected pieces. 




Figure 10: Fermi surface area for neutrons S-^/Sgg^s ^ function of lattice spacing 
(aQy ) below and above phase transition from lasagna (slab) regime to spaghetti (rod) 
regime. 
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turbation in the density region we considered. Taking into account such term would 
raise some degeneracies. The most dramatic change concerns some regions in momen- 
tum space around which unperturbed bands are crossing. The breaking of the spin 
symmetry will entail that such configuration may be turned into band "kissing" (see 
figure El) which means that the group velocity will be strongly reduced around those 
momenta. This means that the resulting mobility would be lower than the one we 
have found. Besides, since the conduction neutron density is rather unsensitive to such 
change, it would not be much affected and therefore the effective neutron mass would 
be slightly larger than the one calculated here. 

6 Conclusion 

The scattering of dripped neutrons by the nuclei in the inner crust leads on a macro- 
scopic scale to a modification of the neutron mass , which can be expressed via a 
well defined mobility scalar )C by = n/JC in which n is the (arbitrary) density 
of such unbound neutrons, and JC is found to be expressible as an integral of the group 
velocity over the corresponding Fermi surface. This effective macro mass should 
not be confused with the effective micro mass, relevant for subnuclear scales, which is 
usually found to be smaller than the ordinary neutron mass. 

Bragg scattering of dripped neutrons is taken into account here by applying Bloch 
type boundary conditions, which are well known in solid state physics but have been 
barely used in this nuclear context. We have computed numerical values of this mobility 
scalar in the bottom layers of the inner crust near the crust-core interface, for simple 
models in the "pasta" layers: equally spaced slab shaped nuclei ( "lasagna" ) and rod like 
nuclei on either a square or an hexagonal 2D lattice ("spaghetti"). The (anisotropic) 
entrainment effect is small at such densities since the system is nearly homogeneous. 
It appears that the mobility scalar tends to be systematically reduced compared to the 
homogeneous expression and the farther from homogeneity, the smaller is the mobility 
scalar. The resulting effective mass 777.^ is found to be larger than the bare neutron 
mass. This results can be interpreted as a macroscopic manifestation of the modifica- 
tions in the shape of the Fermi surface. Notably as one goes from the homogeneous 
outer core to the crust, the spherical neutron Fermi surface gets distorted and even torn 
and pierced. 

The main reason is that the neutron Bragg scattering by crustal nuclei leads to 
the opening of energy band gaps. A specific feature of those exotic phases is that the 
mobility scalar is bounded, JC > 2n^/3m for the "lasagna" phase and JC > n^/3m 
for the "spaghetti" phase due to the fact that the neutrons are still free to move in one 
or two dimensions, respectively. There is no such bound for three dimensional crystals, 
for which smaller mobility scalars (larger effective masses) may be expected. 

From these considerations, we can infer that the mobility scalar JC will increase 
with the density, starting from zero in the outer crust below the neutron drip threshold 
where all neutrons are confined, since = for all states, to its largest possible value 
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in the homogeneous neutron star mantle. At low densities near neutron drip, where 
the conduction neutron density is negligible compared to the total neutron density, 
n <^ J the crystal potential V is quite large but the fraction of the cell occupied by 
the nuclei is small so that dripped neutrons propagate essentially freely. Consequently 
we may expect to get JC — n/m and — m (on the understanding that, as 
discussed in Subsection 13.61 a conduction state is defined so as to have an associated 
group velocity that differs significantly from zero). On the other hand, at very high 
densities where nuclei nearly merge, the potential is much weaker and smoothly varying 
and we shall have n ~ TT-n , so we expect to find ~ m . 

In other words the mean effective mass 772^ is expected to be close to the ordinary 
mass HI at the top (above neutron drip density where neutrons start to leak out of nu- 
clei) and bottom (where nuclei merge into a uniform mixture of neutrons, protons and 
electrons at density about 1/3 or 2/3 the nuclear saturation density ng^t — 0.16 fm 

[25) of the inner crust, but may reach a much larger value (with astrophysically in- 
teresting consequences) in the intermediate layers. This indicates that the unbound 
neutron band effects (that are effectively neglected in the commonly used W-S approx- 
imation) could have important consequences as concerns the equilibrium neutron star 
crust structure and composition. 

The exploratory character of the present work, justifies the simplicity of the single 
particle model we used. In a more realistic, Hartree-Fock calculation, the single particle 
potential and effective micro mass should have to be determined self-consistently. The 
potential may contain momentum dependent terms, as with effective nucleon-nucleon 
interactions of the Skyrme type, which leads to space varying effective neutron masses 
m®{r} which are typically smaller than the bare neutron mass. However the general 
arguments we developed, relying on the shape of the Fermi surface, suggest that the 
qualitative results of our calculations, namely the enhancement of the macroscopic 
effective mass will remain valid in more elaborate single particle schemes. Finally 
let us mention that as we have shown in a recent paper the neutron pairing is not 
expected to qualitatively alter our present conclusions. 
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APPENDIX 

A Numerical results 



A.l Slab shaped nuclei 



nn(fm 3) 


^(MeV) 


aoy(fm) 






m^/m 




0.0735 


31.70 


23.71 


0.9634 


1.0698 


1.0307 


0.9878 


0.0749 


32.10 


23.07 


0.9646 


1.0664 


1.0286 


0.9885 


0.0773 


32.79 


22.23 


0.9666 


1.0605 


1.0251 


0.9896 


0.0792 


33.36 


21.84 


0.9687 


1.0526 


1.0196 


0.9910 



A. 2 Rod like nuclei 

A. 2.1 hexagonal lattice 



^n(fm ^) 


^z(MeV) 


aoy(fm) 


n/nn 




m;^ / m 




0.0581 


27.2461 


27.17 


0.9553 


1.4102 


1.3471 


0.8898 


0.0630 


28.7422 


25.77 


0.9578 


1.2972 


1.2425 


0.9145 


0.0678 


30.1836 


24.62 


0.9606 


1.2225 


1.1744 


0.9322 


0.0716 


31.2812 


23.97 


0.9632 


1.1668 


1.1239 


0.9471 



A. 2. 2 square lattice 



nn(fm ^) 




aoy(fm) 






m^/m 




0.0581 


27.2461 


27.17 


0.9553 


1.4673 


1.4016 


0.8778 


0.0630 


28.7422 


25.77 


0.9578 


1.3231 


1.2673 


0.9085 


0.0678 


30.1797 


24.62 


0.9606 


1.2124 


1.1647 


0.9347 


0.0716 


31.2812 


23.97 


0.9632 


1.1473 


1.1051 


0.9524 
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